Detecting genuine multipartite continuous-variable entanglement 
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We derive necessary conditions in terms of the variances of position and momentum linear com- 
binations for all kinds of separability of a multi-party multi-mode continuous-variable state. Their 
violations can be sufficient for genuine multipartite entanglement, provided the combinations con- 
tain both conjugate variables of all modes. Hence a complete state determination, for example by 
detecting the entire correlation matrix of a Gaussian state, is not needed. 
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I. INTRODUCTION 



Quantum entanglement shared by two parties en- 
hances their capability to communicate. In principle, 
it allows them to convey quantum information reliably 
through a classical channel (quantum teleportation 0), 
to double the amount of classical information transmit- 
table through a classical channel (quantum dense coding 
0), or to prepare at a distance states from nonorthogo- 
nal bases for secure communication (quantum key distri- 
bution H H> Hi)- These entanglement-assisted communi- 
cation schemes are extendible to an arbitrary number of 
parties sharing multipartite entanglement. For instance, 
a sender may transfer quantum information through clas- 
sical channels to several receivers as reliably as allowed 
by optimal cloning (telecloning Q), or the parties may 
share quantum (or classical) information retrievable only 
when all parties cooperate (quantum secret sharing 0). 
A more recent proposal that exploits the multi-party 
quantum correlations of multipartite entangled states is 
the so-called Byzantine agreement protocol In gen- 
eral, the presence of entanglement is verified through the 
success of a quantum protocol that would fail otherwise 
(e.g., quantum teleportation). Such an operational crite- 
rion is only sufficient for entanglement and failure does 
not necessarily rule out its presence. In an experimen- 
tal realization, however, before running through an en- 
tire entanglement-based protocol, it is desirable first to 
confirm that the generation of sophisticated multi-party 
entangled states has succeeded. The aim of this paper is 
to provide a simple but unambiguous experimental test 
to check for a particular kind of genuinely multipartite 
entangled states, namely those described by continuous 
variables (cv) and produced with squeezed light and lin- 
ear optics. 

Work in the direction of generating tripartite cv en- 
tanglement has been carried out already by sending one 
half of a two-mode two-party entangled state through 
an extra beam splitter with a coherent state or a vacuum 
state at its second input port @, . The resulting three- 
mode state was a side product of the Bell measurement 
for the teleportation of coherent states using a preshared 
symmetric two- mode squeezed state H . Its tripartite en- 



tanglement was not further investigated in that experi- 
ment. In another experiment, reported recently [fUj , the 
two-mode state was asymmetric, corresponding to two 
independently squeezed states combined at an asymmet- 
ric 1 : 2 beam splitter. The output three-mode state after 
an additional symmetric beam splitter was then similar 
to the states proposed in Ref. Q. 

Quantum communication, or more general quantum 
information with cv has attracted a lot of interest due 
to the relative simplicity and high efficiency in the gen- 
eration, manipulation, and detection of optical cv states. 
Although recent results suggest that these assets of Gaus- 
sian cv operations (phase shifting, beam splitting, homo- 
dyne detections, phase-space displacements, squeezing) 
are not extendible to more advanced quantum protocols 
such as entanglement distillation [l4| , the simple 

and efficient cv approach still seems promising for many 
tasks and might be suitable for others too when combined 
with discrete- variable (dv) strategies. On the other hand, 
potential linear-optics implementations of quantum pro- 
tocols solely based on dv utilizing single photons are re- 
stricted by No-Go results such as the impossibility of a 
complete distinction between the four Bell states . In 
order to perform such a Bell measurement near perfectly 
with linear optics, one has to employ complicated entan- 
gled states of many auxiliary photons [Q. In contrast, 
a Bell and also a GHZ state analyzer can be easily con- 
structed in the cv setting using only beam splitters and 
homodyne detectors ftnj ]l8|, [L9| . 

How may one now verify experimentally the presence 
of entanglement without implementing a full quantum 
protocol? We are here particularly concerned about the 
experimental verification of genuinely multipartite en- 
tangled states where none of the parties is separable 
from the rest (in terms of the separability properties of 
the total density matrix). In general, theoretical tests 
might be as well applicable to the experimental verifica- 
tion. For instance, the violation of inequalities imposed 
by local realism confirms the presence of entanglement. 
Proving genuine multipartite entanglement, however, re- 
quires stronger violations 1 20 than those determined by 
the commonly used Mermin-Klyshko TV-party inequal- 
ities [pll Moreover, in any case, violations of Bell- 
type inequalities using Gaussian cv entangled states with 
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always positive Wigner functions must rely on observ- 
ables other than the quadratures (i.e., position and mo- 
mentum). Photon number parity may serve as an ap- 
propriate dichotomic variable to reveal the nonlocality 
of the cv entangled states pal . This applies to the two- 
party two-mode EPR-like p3] and to the iV-party TV- 
mode GHZ-like cv states [B4L Such an approach, how- 
ever, is not very feasible due to its need for detectors 
resolving large photon number. 

The negative partial transpose (npt) criterion |2[| is 
sufficient and necessary for the bipartite inseparability 
of 2 x 2-dimensional, 2 x 3-dimensional fiij, and 1 x JV 
mode Gaussian states 
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2q |. A complete experimental 
determination of the state in question would also enable 
an npt check. In general, any theoretical test is applica- 
ble when the experimentalist has full information about 
the quantum state after measurements on an ensemble 
of identically prepared states (e.g. by quantum tomogra- 
phy ]29|l). Such a direct verification of entanglement via a 
complete state measurement is in general very demand- 
ing to the experimentalist, in particular when the state 
to be determined is a potentially multi-party entangled 
multi-mode state. 



II. GAUSSIAN STATES 

The multi-party entanglement criteria that we will de- 
rive here do not rely on the assumption of Gaussian 
states. However, the states commonly produced in the 
laboratory are indeed Gaussian and the theoretical clas- 
sification of different types of multipartite entanglement 
becomes simpler for Gaussian states |50| . 

Since the entanglement properties of a multi-mode 
multi-party state are invariant under local phase-space 
displacements, the multi-mode states may have zero 
mean and their Wigner function is of the form 



W(£) 



1 
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exp 



(i) 



with the 2,/V-dimcnsional vector £ having the quadrature 
pairs of all N modes as its components, 

£ = {xi,pi,x 2 ,p2, — ,xn,Pn) , (2) 

and with the 2N x 2N correlation matrix V"W having as 
its elements the second moments symmetrized according 
to the Weyl correspondence 

Tr[p + A&A&)/2] - ((Mj + &ii)/2) 



(JV) 



(3) 



where A£, = — (£,•) = £j for zero mean values. Note 
that the correlation matrix of any physical state must be 



real, symmetric, positive, and must obey the commuta- 
tion relation (g7|, , 



A 



kl , 



k,l= l,2,3,...,2iV, 



(4) 



with the 2A x 2N matrix A having the 2x2 matrix J 
as diagonal entry for each quadrature pair, for example 
for N = 2, ' 



A = 



J 
J 



J = 



1 

-1 



(5) 



A direct consequence of this commutation relation and 
the non-negativity of the density operator p is the N- 
mode uncertainty relation p7L Pq| , 



F W_1 A >0 
4 ~ 



(6) 



Note that this condition is equivalent to +iA/4 > 
by complex conjugation. As for the direct verification 
of entanglement via a complete state measurement, for 
Gaussian cv states, the complete measurement of an N- 
party A-mode quantum state is accomplished by deter- 
mining the 2N x 2N second-moment correlation matrix. 
This corresponds to N(l + 2N) independent entries tak- 
ing into account the symmetry of the correlation matrix. 
Kim et al. |52j recently demonstrated how to determine 
all these entries in the two-party two-mode case using 
beam splitters and homodyne detectors. Joint homo- 
dyne detections of the two modes yield the intermode 
correlations such as (X1X2) — {xi){x2), (x\P2) — {xx)(p2), 
etc. Determining the local intramode correlations such 
as {xipi +pxXx)/2 — (xi)(pi) is more subtle and requires 
additional beam splitters and homodyne detections (or, 
alternatively, heterodyne detections). Once the 4x4 
two-mode correlation matrix is known, the npt criterion 
can be applied as a sufficient and necessary condition for 
bipartite Gaussian two-mode inseparability (where npt 
corresponds to a sign change of all the momentum vari- 
ables with positions unchanged pTJ]). In fact, the en- 
tanglement can also be quantified for a given correlation 
matrix |}2| [33]. For three-party three-mode Gaussian 
states, one may pursue a similar strategy. After measur- 
ing the 21 independent entries of the correlation matrix 
(for example, by extending Kim et al.'s scheme (3^] to 
the three-mode case), the sufficient and necessary crite- 
ria by Giedke et al. |30| can be applied. Let us examine 
the separability properties of (in particular, three-party 
three-mode) Gaussian states in more detail. 



III. SEPARABILITY PROPERTIES OF 
GAUSSIAN STATES 

The criteria by Giedke et al. determine to which of 
five possible classes of fully and partially separable, and 
fully inseparable states a three-party three-mode Gaus- 
sian state belongs. Hence genuine tripartite entangle- 
ment if present can be unambiguously identified. The 
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classification is mainly based on the npt criterion for 
cv states. Transposition is a positive map that corre- 
sponds in phase space to a sign change of all momentum 
variables, £ T -> T£ T = (xi, -pi, x 2 , -P2, —, xn, ~Pn) T 
p7[ . In terms of the correlation matrix, we have then 
V^™) — * rV^r. Since transposition is not a completely 
positive map, its partial application to a subsystem only 
may yield an unphysical state when the subsystem was 
entangled to other subsystems. Expressing partial trans- 
position of a bipartite Gaussian system by T a = T ® II 
(where A® B means the block-diagonal matrix with the 
matrices A and B as diagonal 'entries', and A and B 
are respectively 2N x 2N and 2M x 2M square matrices 
applicable to N modes at a's side and M modes at 6's 
side), the condition that the partially transposed Gaus- 
sian state described by T a V^ N+M ^T a is unphysical [see 
Eq. d)], r a V<- N+M) r a £ | A, is sufficient for the insep- 
arability between a and b [^7[ |28). For Gaussian states 
with N = 1 and arbitrary M , this condition is sufficient 
and necessary |2^] . The simplest example where the con- 
dition is no longer necessary for inseparability involves 
two modes at each side, N = M = 2. In that case, states 
with positive partial transpose (bound entangled Gaus- 
sian states) exist p8fl . For the general bipartite N x M 
case, there is also a sufficient and necessary condition: 
the correlation matrix V^ N+M ' corresponds to a sepa- 
rable state iff a pair of correlation matrices Va N ^ and 



V h {M > exists such that V^ N+M '> > V} 



(AT) 



■V b (M) . Since it 
is in general hard to find such a pair of correlation ma- 
trices Vo and V^ M ^ for a separable state or to prove 
the non-existence of such a pair for an inseparable state, 
this criterion in not very practical. A more practical 
solution was provided in Ref. (34). The operational cri- 
teria there, computable and testable via a finite number 
of iterations, are entirely independent of the npt crite- 
rion. They rely on a nonlinear map between the corre- 
lation matrices rather than a linear one such as the par- 
tial transposition, and in contrast to the npt criterion, 
they witness also the inseparability of bound entangled 
states. Thus, the separability problem for bipartite Gaus- 
sian states with arbitrarily many modes at each side is 
completely solved. For three-party three-mode Gaussian 
states, the only partially separable forms are those with 
a bipartite splitting of 1 x 2 modes. Hence already the 
npt criterion is sufficient and necessary. 

The classification of tripartite three-mode Gaussian 
states l30f, 



class 1 : 
class 2 : 
class 3 : 
class 4 or 5 : 



v, 



(3) 



> 



(3) 



> 4 A ^3 3) > i A- 



- (3) 

is solely based on the npt criterion, where = 
TjV^Tj denotes the partial transposition with respect 
to one mode j. In classes 2 and 3, any permutation 
of modes (k, m, n) must be considered. Class 1 cor- 
responds to the fully inseparable states. Class 5 shall 
contain the fully separable states. For the fully separa- 
ble Gaussian states if described by V^, one-mode cor- 
relation matrices V x , V% 1 \ and V3 exist such that 



(i) 



V, 



(i) 



VS 1 ' . In general, fully separable 



V< 3 > > vi 

quantum states can be written as a mixture of tripar- 
tite product states, J2iViPi,i ® Pi,2 ® Pi, 3- In class 2, 
we have the one-mode biseparable states, where only one 
particular mode is separable from the remaining pair of 
modes. This means in the Gaussian case that only for 

one particular mode k, > (BVmn with some two- 

(2) 

mode correlation matrix Vmn and one-mode correlation 
matrix Vu . In general, such a state can be written as 
J2i Vi Pi,k ® Pi,mn for one mode k. Class 3 contains those 
states where two but not three bipartite splittings are 
possible, i.e., two different modes k and m are separable 
from the remaining pair of modes (two-mode bisepara- 
ble states). The states of class 4 (three-mode biseparable 
states) can be written as a mixture of products between 
any mode 1, 2, or 3 and the remaining pair of modes, 
but not as a mixture of three-mode product states. Ob- 
viously classes 4 and 5 are not distinguishable via the 
npt criterion. An additional criterion for this distinction 
of class 4 and 5 Gaussian states is given in Ref. p(J , de- 

(rTV(i) 



ciding whether one-mode correlation matrices v£ 1 , VC 



(7) 



and Vf> exist such that > 1 ' Vf>. For 
the identification of genuinely tripartite entangled Gaus- 
sian states, only class 1 has to be distinguished from the 
rest. Hence the npt criterion alone suffices. 

What about more than three parties and modes? Even 
for only four parties and modes, the separability issue 
becomes more subtle. The one-mode bipartite splittings, 
J2i Vi Pi,Mm <8 Pi, n > can be tested and possibly ruled out 
via the npt criterion with respect to any mode n. In the 
Gaussian language, if Vn^ 1 ^ 4 A for any n, the state can- 
not be written in the above form. Since we consider here 
the bipartite splitting of 1 x 3 modes, the npt condition 
is sufficient and necessary for Gaussian states. However, 
also a state of the form r\i pi,u®Pi,mn leads to negative 
partial transpose with respect to any of the four modes 
when the two pairs (k, I) and (m, n) are each entangled. 
Thus, npt with respect to any individual mode is nec- 
essary but not sufficient for genuine four-party entangle- 
ment. One has to consider also the partial transposition 
with respect to any pair of modes. For this 2x2 mode 
case, however, we know that entangled Gaussian states 
with positive partial transpose exist |2^]. But the npt 
criterion is still sufficient for the inseparability between 
any two pairs. As for a sufficient and necessary condition, 
one can use those from Ref. p4| . In any case, in order 
to confirm genuine four-party or even TV-party entangle- 
ment, one has to rule out any possible partially separable 
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form. In principle, this can be done by considering all 
possible bipartite splittings (or groupings) and applying 
either the npt criterion or the stronger operational crite- 
ria from Ref. Hj. Although a full theoretical character- 
ization including criteria for entanglement classification 
has not been considered yet for more than three parties 
and modes, the presence of genuine multipartite entan- 
glement can be confirmed, once the complete 2N x 2N 
correlation matrix is given. 

Rather than detecting all the entries of the correlation 
matrix we are aiming here at a simple check based on 
only a few measurements, preferably efficient homodyne 
detections. Even for larger numbers of parties, this check 
should remain simple. Though it may not yield full infor- 
mation (e.g., the complete correlation matrix) about the 
quantum state of interest, it should still unambiguously 
verify the presence of genuine multipartite entanglement. 
This check may prove the presence of entanglement indi- 
rectly through measurements after transforming the rele- 
vant state first into an appropriate form via linear optics. 



IV. DETECTING ENTANGLEMENT: 
BIPARTITE CASE 

In the two-party two-mode case, the necessary separa- 
bility condition for any cv state p5| 

<[A(ir -i 2 )] 2 ) + <[A(p! +p 2 )] 2 ) > 2\{%p])\ , (8) 

can be tested, for example, with a single beam splitter. 
The position and momentum variables x\ and pi (units- 
free with H = h, [%i,Pk] = correspond to the 
quadratures of two electromagnetic modes, i.e., the real 
and imaginary parts of the annihilation operators of the 
two modes: ai = x\ + ipi. The beam splitter provides 
the suitable quadrature combinations for the positions 
and momenta simultaneously detectable at the two out- 
put ports. Without beam splitter, just by measuring 
first both positions and subtracting them electronically, 
and in a second step detect ing both momenta and com- 
bining these electronically [p6| , a more direct test of the 
two-party condition is possible. However, instead of a 
simultaneous detection of the relevant combinations, it 
requires switching the two local oscillator phases from 
position to momentum measurements. For an ensemble 
of identically prepared states, this sequence of detections 
would still enable the application of the two-party con- 
dition. Note that the violation of Eq. (||) is only suffi- 
cient for inseparability, i.e., there are (even Gaussian) cv 
entangled states that satisfy Eq. (||). Any Gaussian cv 
state, however, can be transformed via local operations 
into a standard form and the presence of entanglement 
would then always yield a violation (3^] (alternatively, 
one may modify the inequality and leave the Gaussian 
state unchanged to obtain a sufficient and necessary con- 
dition p^]). The point is that the entanglement of states 
already in this standard form (such as two-mode squeezed 



states) can, in principle, always (for any nonzero squeez- 
ing) be verified experimentally by checking Eq. (g). A 
full determination of the correlation matrix, including 
elements such as (xipz) — {xi)(p 2 ) which do not appear 
in the expressions of Eq. (B|), is not required. Measur- 
ing also these elements may confirm that the state is in 
standard form (when they are zero) and hence render the 
condition Eq. (g) sufficient and necessary for separabil- 
ity. In any case, it would also enable quantification of 
the entanglement |3^, 33 . 

The combinations in condition Eq. (^) are exactly 
those detected in a cv Bell measurement of modes 1 and 
2 | p7[ . Thus, the verification of non-maximum two-mode 
cv entanglement may rely on measurements of observ- 
ables that are detected for the projection onto the max- 
imally entangled cv basis of two modes. Now we investi- 
gate the TV-party N-mode case in that respect. 



V. THE CV GHZ BASIS 

Let us introduce the maximally entangled states 

1 f°° 

|*(u,ui,u 2 ,-,«Ar-i)> = —;= \ dxe 2lvx 

V 7 *" J-oo 

x | a:) ® | a; — Wi) ® \x — u\ — u 2 ) 

<g> • <g> |x - ttl - U 2 MJV-l) ■ (9) 

Since f^L. \x){x\ — lL and (x\x r ) — 8(x — x'), they form 
a complete, 

dv du±du2- • • dupf-i (10) 

X\V(V,U1,U2, ...,U N - 1 ))(^(v,Ui,U 2 , ...,UjV-l)| = ti-® N , 

and orthogonal, 

(*(U,U1,U2, ..^UjV-l)!*^,"!,^! -,u' N -l)) (11) 

= 5(v- v')S(ui - u[)S(u 2 - u'2) ■ ■ -S(u N -i - Ujy_i), 

set of basis states for N modes. In a "cv GHZ state 
analyzer" , determining the quantities v = p\ + p 2 + 
•■• + p N , ui = xi — x 2 , u 2 = x 2 - x 3 ,..., and 
uat_i = xn-i — xn means projecting onto the basis 
{|^(u, tti, U2, kjv-i))}- This can be accomplish with 
a sequence of beam splitters and homodyne detections 
JlSj, |l9| |. Inferring from the two-party case, we may con- 
jecture that the N quadrature combinations given by 
v, tti, U2, ...,itjV-i provide a sufficient set of observables 
for the verification of (possibly genuine) iV-party entan- 
glement. Just as for two parties, the variances of these 
quantities could then also be determined by combining 
the results of direct x and p measurements electronically. 
It was shown in Ref. [[l8| [l9| that conditions for genuine 
multipartite entanglement can be derived based on the 
above N combinations and additional assumptions such 
as the purity and the total symmetry of the state in ques- 
tion. Later we derive a set of A — 1 conditions for those N 
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combinations sufficient for the presence of genuine mul- 
tipartite entanglement. This set is well suited for the 
experimental confirmation of the genuine multi-party en- 
tanglement of cv GHZ-type states. No extra assumptions 
about the state are needed in order to close the loophole 
of partial separability. First, we discuss now what the 
structure of simple experimental criteria for multipartite 
cv entanglement might be. 



VI. DETECTING ENTANGLEMENT: 
TRIPARTITE CASE 

Let us consider three parties and modes. The goal is 
to extend the simple two-party two-mode entanglement 
check to a simple test for genuine three-party three-mode 
entanglement. The criteria are to be expressed in terms 
of the variances of quadrature linear combinations for the 
modes involved. Defining 

u = h x xi + h 2 x 2 + h 3 x 3 , v = gip! + g 2 p 2 + 93P3 , (12) 

a fairly general ansatz is 

((Au) 2 ) p + ((Av) 2 ) p > f(hi,ha,h3, 9i, 52, 33) , (13) 

as a potential necessary condition for an at least partially 
separable state. The position and momentum variables 
xi and pi are the quadratures of the three electromagnetic 
modes. The hi and gi are arbitrary real parameters. We 
will prove the following statements for (at least partially) 
separable states, 

P = E ^ Pi,n <8 pi,3 

i 

- f(hi,gi) = {\h a g 3 \ + \h 1 g 1 + h 2 gi\)/2, (14) 

p = e Vi Pi.ia ® Pt,a 

i 

-» f{h h gi) = {\h 2 g 2 \ + \h l9l + h 3 g 3 \)/2, (15) 

p = e^ Pi, 23 <?> pi,i 

i 

-> f(h U gi) = (\h l9l \ + \h 2 g 2 + h 3 g 3 \)/2. (16) 

Here, for instance, Pi,i 2 (8 Pi,3 indicates that the three- 
party density operator is a mixture of states i where par- 
ties (modes) 1 and 2 may be entangled or not, but party 
3 is not entangled with the rest. Hence also the fully 
separable state is included in the above statements. In 
fact, for the fully separable state, we have 

P = E ^ Pi,l ® Pi,2 ® Pi,3 
i 

f{h h gi) = (\h m \ + \h 2 g 2 \ + \h 3 g 3 \) /2,(17) 

which is always greater or equal than any of the bound- 
aries in Eq. Q, Eq. @, or Eq. @. For the proof, let 
us assume that the relevant state can be written as 



For the combinations in Eq. (|T^), we find 

((Au) 2 ) p + {(Av) 2 ) p 
= ({u% + {v\)-(u) 2 p -{v)l 

i 

= E 7 ^ \hl{xl)i + Ki{^m)i+K{ & V}i 

i 

+9k{Pk)i + 9 2 m(p 2 m)i + 9l{pt)i 

~~^2 (^h k h m (x k x m ^ i -\- h k h n ( K x k x n ) i -f- h rn h n (^x rn x n ^ 

+2(gk9m(pkPm)i + 9k9niPkPn)i + 9m9n{PmPn) ij 
-(u) 2 p -(v) 2 p 

= ^Vi [^((Ai fc ) 2 ) 4 + / l 2 m ((Ai„ i ) 2 ) 4 + ^((Ai n ) 2 ) 

i 

+g 2 ((Ap k ) 2 ) l +g 2 m ((Ap m ) 2 ) l +g 2 n ((Ap n ) 2 ) 

+2h k h n (^(x k x n )i - (x k )i(x n ) 
~\~2,h rn h n ^{x m x n ^>i 
+2g k g in ((p kPm/i 

+2gk9n({PkPm)i - (pk)i{p m )i 

+2g m g n (<P 

mPn/i (pm)i(Pn) 
+ E'? 4 &)i ~ E % (u)i 



E 



Vi w 



(19) 



Vi Pi, km yy Pi,n 



(18) 



where (■■■){ means the average in the state pi ykm ® Pi, n - 
Note that in the derivation so far we have not used the 
particular form in Eq. ( |l8|) yet. Exploiting this form of 
the state, we obtain (x k x n )i = {x k )i(x n ) u (x m x n )i = 
(x m )i(x n )i, and similarly for the terms involving p. Be- 
cause modes k and m may be entangled in the states i, 
we cannot replace (x k x m )i by (x k )i{x m )i, etc. By apply- 
ing the Cauchy-Schwarz inequality as in the two-party 
derivation of Rcf. ||, Ei-WL^ (EiW)il) , we 
see that the last two lines in Eq. ([19]) are bounded below 
by zero. Hence in order to prove (jAu) 2 ) p + ((Av) 2 ) p > 
(\h n g n \ + \h k g k + h m g m \)/2, it remains to be shown that 
for any i [recall that the mixture in Eq. jl^ ) is a convex 
sum with E 4 Vi = 1] 1 

h 2 {(Ax k ) 2 ) t + h 2 n {(Ax m ) 2 ) t + h 2 n ((Ax n ) 2 )i 
+g 2 ((Ap k ) 2 ) l +g 2 m ((Ap m ) 2 ) l +gl((Ap n )% 

+2g k g 

((p kPm/i (pk)i(p 
> (\Kgn\ + \h k g k + h m g m \)/2 . (20) 
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By rewriting the left-hand-side of Eq. (|2' 
variances only, indeed we find 



in terms of 



+ ([A(h k x k + h m x m )] 2 ) l + ([A(g k p k + 

!JmPm ) J )i 

> \([Kx n , 9nPn])\ + \ ([hk 

= (\h n g n \ + \hk9k + h m g m \)/2 , (21) 

using the sum uncertainty relation ((AA) 2 ) + ((AB) 2 } > 
\{\A, &\)\ and [xi,pj] = iSij/2. Hence the statements in 
Eq. (UJ) are proven when we consider the correspond- 
ing permutations of (fc, m, n) = (1, 2, 3). The inequalities 
Eq. (|l|) with Eq. @), Eq. ((TsJ) , and Eq. (|l|) represent 
necessary conditions for all kinds of (partial) separabil- 
ity in a tripartite three-mode state. One may then prove 
the presence of genuine tripartite entanglement through 
violations of these inequalities, thus ruling out any (par- 
tially) separable form. Whether there are really three 
different conditions required for the verification depends 
on the choice of the coefficients hi and g\ in the linear 
combinations. For a particular choice, some of the con- 
ditions may coincide. For example, consider hi = fli = 1 
and <?2 = 33 = — fa = —^3 = 1 /y/2 in Eq. (12). In this 
case, the boundaries in Eq. ( |l4| ) and Eq. (15) become 
identical, f[hi,g{) = 1/2. The boundary of Eq. (16) is 
even larger, f(hi,gi) = L equivalent to that for a fully 
separable state in Eq. ( |l7| ) . Hence the violation of a sin- 
gle condition, 

({A[£ 1 -(x 2 +£ 3 )/V2]} 2 ) p 

+ ({A[p 1 + (p 2 +p 3 )/V2]} 2 ) p > 1/2, (22) 

is already sufficient for genuine tripartite entanglement. 
These particular combinations are not only significant 
for the reason that they yield nonzero boundaries for all 
kinds of separable states. Moreover, their commutator 
vanishes, 

[xi-(i2+X3)/%/2,pi + (p 2 +p 3 )/V2] = 0, (23) 



allowing for arbitrarily good violations of Eq. (|22|) and, 
in principle, the existence of a simultaneous eigenstate of 
these two combinations. Such a state corresponds to the 
three-mode state obtainable by splitting one half of an 
infinitely squeezed two-mode squeezed (EPR) state at a 
50:50 beam splitter. The EPR correlations, x\ — x 2 — > 
and pi + p2 — * 0, are then transformed into the three- 
mode correlations x\ — (x' 2 +x' 3 )/ V2 — > and pi + (p' 2 + 
p' 3 )/y/2 -► 0. Let us turn to an arbitrary number of 
parties (modes) now. 



the following. First, measure both quadratures x and p 
of all modes involved and combine them in an appropri- 
ate linear combination. The variances of these combina- 
tions may then yield violations of conditions necessary 
for partial separability. Appropriate combinations are 
those where the total variances for all partially separable 
states have nonzero lower bounds and where the commu- 
tators of the combinations vanish. As for the derivation 
of the corresponding entanglement criteria, we employ 
the following steps. 

1. Select a distinct pair of modes (m,n). 

2. Choose appropriate linear combinations of the 
quadratures in order to rule out all possible separable 
splittings between this pair of modes in the convex sum 
of the total density operator. 

3. Consider different pairs (m, n) to negate all partial 
separabilities; if necessary add further conditions involv- 
ing other linear combinations. 

Below it will become clear that step 2 can be performed 
simply by using the appropriate bipartite combinations, 
x m - x n and p m + p„, i.e., by taking all hi = gi = 
except h m = g m = 1 and h n = — g n = — 1 in the general 
combinations 



u — hiXi + h 2 X2 + ■ ■ ■ + h^XN 
v = gipi + g 2 p2 H h 9nPn ■ 



(24) 



The boundaries of the total variance conditions are then 
identical for any pair (to, n) separable in the convex sum, 
namely f(hi,gi) = 1 in 



((Au) 2 ) p + ((Av) 2 ) p > f(hi,h 2 ,...,hN,9i,g2,-,9N) 



(25) 



However, in general, one obtains better multi-party con- 
ditions when linear combinations for the quadratures of 
more than only two modes are used. Through such multi- 
mode combinations the potential multi-mode correlations 
are taken into account. Before giving an example, let us 
first derive the general A^-party bounds in the condition 
Eq. (|2^) . For any partially separable form, the total den- 
sity operator can be written as 



VII. DETECTING ENTANGLEMENT: 
MULTIPARTITE CASE 



Pi,k 3 



(26) 



Inferring from the discussion of the previous section, 
the recipe for verifying the genuine multipartite entan- 
glement between arbitrarily many parties and modes is 



with a distinct pair of "separable modes" (to, n) and the 
other modes k r ^ k s . For the combinations in Eq. 
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we find now 

((Au) 2 ) p + ((Av) 2 ) p 

= E^ «* 



E 7 ^ \h 2 n {x 2 m ) l + h 2 n {x 2 n/l 



N-2 



J2 h l 

3=1 

JV-2 

+5m(Pm)i +9l(pl)i + E 
3=1 

X! )' + flMfe,/ (Pfc 3 

3W=1 
JV-2 

3=1 

Pm)i H" 9kj9n {PkjPn) 
-t~2 (^h m h n (x m X n ) ^ -(- Q-rnQniPmPn) 

-{u)l-{v? P 



N-2 



Pk H 



= {^„((A£ m ) 2 > 2 + ^((Ai 

^ 1 ^((Ap 



n v ) i 



9m(( A Pm/ H 

N-2 

E (^,(( A -^) 2 )«+5M( A ^) 2 



3 = 1 



)i( x k r ,) 

+9k r 9k r , ({Pk r Pk r ,)i - (Pk r }i{Pk r ,)i) 

+9k s gk B , ((pk 3 Pk s ,)i - (pk a )i(pk a ,) 

+2 E [^fcr^mU 
r 

+9k r g m \{Pk r P m )i ~ (Pk r )i(Pm) 

%k s &n) i {%ks)i (%ri)i 
s 

+9k a gn((pk,Pn)i ~ (Pk 3 ) l (Pn)?j } 

+E 7?i ^ 2 ~ ^E^ 

+e^^ ? -(e^^m • 



(27) 



For the last equality, we exploited Eq. (£6|), namely 
that modes k r through m are separable from modes k s 
through n in the convex sum of the total density oper- 
ator. Similar to the three-party case, we can now apply 



the Cauchy-Schwarz inequality to the last two lines of 
Eq. ( p7j ) and express the remaining terms by variances 
only. This leads for any i to 



+ 



A\h m x rn + E h k r x kr 

r 

+ (^ A (g m p m + E 9k r Pk r 

r 

/\{h n x n + E h ks x ks 

S 

A^gnPn+^gkJk^ 

s 

hmXm + hKXkr. , 9mPm + 9k r Pk T 
r r 

h n X n + E h k B X ks , gnPn + E 9k s Pk 3 



> 



(28) 

using again the sum uncertainty relation ((AA) 2 ) + 
((Av3) 2 ) > Thus, by evaluating the commuta- 

tors with = iSij/2, we obtain for the total variance 



<(A«) 2 )p 
1 



<(A£) 2 } P 



> 



hrn9n 



h kr g kr + h n g n + E h 



k 3 gk s 



(29) 



Any additional splitting of the parties in the states i, 

P = J2t Vi Pi,k r -m ® • • • ® Pi,fe r / ® Pi,k 3 -n ® • • • ® , 

would in general make the bound larger, eventually yield- 
ing the bound for the fully separable state, J2j l^jQj]/^ 
(j = l..JV). 

As mentioned previously, the well-known bipartite 
combinations applied to modes (m, n), x m — x n and 
Pm + Pn, mean all hi = gi = except /i m = p m = 1 
and h n — —g n = —1 in Eq. and hence ((Au) 2 ) p + 
<(A«) 2 ) P > 1. 

As for a simple example, we may extend that from the 
previous section to TV modes and set hi = g\ — 1 and 
.92 = .93 = • • • = .9iV = -/i>2 = ~h% = ■ ■ ■ = -h N = 
l/y/N - 1. Without loss of generality, we choose m = 1 
and obtain for a state of the form Eq. (EF 



((A 



it 



1 

> - 
- 2 



1 - 



((A£) 2 ) P 

M r 



N -1 



1 + M s 



N - 1 



(30) 



where M r is the number of modes potentially entangled 
with mode m = 1 in the convex sum and M s is the num- 
ber of modes potentially entangled with mode n in the 
convex sum. Apart from the fully inseparable case M r — 
N — 1, the boundary in Eq. ([30]) is always greater than 
zero allowing for an ultimate nonzero bound for all kinds 
of partial separability. Since [u, v] = 0, genuine iV-party 
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entanglement can be verified when ((Am) 2 ) p + ((At)) 2 ) p is 
sufficiently close to zero. The ultimate (smallest) bound 
is given by the state with the maximum number of modes 
M r inseparable from mode m = 1 in the convex sum, 
M r = TV — 2, and hence M s = 0. This bound is then 
1/(TV — 1). If none of the modes is inseparable from mode 
m = 1, M r = and M s = TV — 2, the boundary becomes 
simply that of a fully separable state, namely one. Thus, 
again the violation of a single condition, 



tion matrix, 



y W = _ 
4 



/ a 



a 



d 
c c 
d d b 



a 



(32) 



where 



{{A[xi - (x 2 + x 3 + ■■ 

+{{A\pi + (p 2 +P3 + 
>l/(N- 



x N )/VW=l]f 



■pn)/Vn=T]} 



1) 



(31) 



is sufficient for genuine TV-partite entanglement. As an 
example for the violation of the ultimate bound for gen- 
uine TV-party entanglement, consider the TV- mode state 
that emerges after symmetrically splitting one half of an 
infinitely squeezed two-mode squeezed state by TV — 2 
beam splitters. The output state is a simultaneous eigen- 
state of u and v. In this case, the EPR correlations, 
x\ — x 2 —> and p\ + p 2 —>■ 0, are transformed into the 
TV- mode correlations x\- 



■_+£' N )/VN^l - 
0. A more sym- 



and pi + (p' 2 +p' 3 + ■■■ +p' N )/^N-l 
metric example is where both halves of an EPR state 
are symmetrically split at beam splitters. The appropri- 
ate combinations are then for instance for four modes, 
x[ + x' 2 — x' 3 — x' 4 and p[ + p' 2 + p' 3 + p' 4 , also having 
zero commutator [u, v] = 0. As a further example, we 
will now discuss the cv GHZ-type states with quadrature 
correlations analogous to those of dv GHZ states. 



VIII. EXAMPLE: CV GHZ-TYPE STATES 



We consider a family of genuinely TV-party entangled 
states. The members of this family are those states that 
emerge from a particular sequence of TV — 1 phase-free 
beam splitters ("TV-splitter") with TV squeezed state in- 
puts Q. By choosing the squeezing direction of one 
distinct input mode orthogonal to that of the remaining 
input modes (mode 1 squeezed in p and the other modes 
squeezed in x, as shown in Figs. @ and | for TV = 3) 
and the degree of squeezing r\ of mode 1 potentially dif- 
ferent from that of the other modes (which are equally 
squeezed by r 2 ) |39(1 , the output states have the following 
properties jl8], |1 Qf| - They are pure TV-mode states, to- 
tally symmetric under interchange of modes, and retain 
the Gaussian character of the input states. Hence they 
are entirely described by their second-moment correla- 



b = — 



d = 





TV- 1 


TV 


TV 


le- 2 - - 


TV- 1 


TV 


" TV 




- e" 2r2 


TV 




-(e- 2ri 


- e +2r2 


TV 





-2r 2 



3 +2r 2 



(33) 



For squeezed vacuum inputs, the multi-mode output 
states have zero mean and their Wigner function is of the 
form Eq. (0). The particularly simple form of the corre- 
lation matrix in Eq. ( |32| ) is, in addition to the general 
correlation matrix properties, symmetric with respect to 
all modes and contains no intermode or intramode x- 
p correlations (hence only four parameters a, 6, c, and d 
are needed to determine the matrix) . However, the states 
of this form are in general biased with respect to x and 
p (a 7^ b). Only for a particular relation between the 
squeezing values (ri,r 2 ) 18, [l9|, 

(TV-1) (34) 



= ±2ri 



x sinh 2?'2 



(TV-1) 2 sinh 2 2r 2 



± 1 



the states are unbiased (all diagonal entries of the corre- 
lation matrix equal), thus having minimum energy at a 
given degree of entanglement or, in other words, max- 
imum entanglement for a given mean photon number 
p9| . The other TV-mode states of the family can be con- 
verted into the minimum-energy state via local squeez- 
ing operations fl|, ||. Only for TV = 2, we obtain 
r = r\ = r 2 . In that case, the matrix V^ N ' reduces to 
that of a two-mode squeezed state which is the maxi- 
mally entangled state of two modes at a given mean en- 
ergy with the correlation matrix entries a = b = cosh2r 
and c = sinh2r = — d. For general TV, the first squeezer 
with r\ and the TV — 1 remaining squeezers with r 2 
have different squeezing. In the limit of large squeezing 

(35) 



(sinh2r2 ~ e + 2 /2), we obtain approximately pi 



(TV-1) 



3 +2r 2 



For the whole family of TV-party TV-mode states with the 
correlation matrix in Eq. ([52]), the quadrature combina- 
tions relevant for detecting genuine multi-party entangle- 
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FIG. 1: Verification of genuine tripartite cv entanglement. 
x measurements: directly detecting the x quadratures of all 
three modes and electronically combining them in an ap- 
propriate way. The three-mode tripartite entangled state of 
modes 1,2, and 3 in this figure is produced with three squeez- 
ers and two beam splitters (the star denotes a 1 : 2 BS). 



merit are [0 || |jJJ 



([A(x m - x n )] 2 ) = e- 2r V2 , 



A \p m +Pn+g [N) J2 Pi 

[2 + (jV-2)gW] 2 c _ 2ri + (, 9 W-l) 2 (jV-2) ^ +2r2 



AN 



2N 



(36) 



The total variances are then optimized (minimized) for 



AN) 

i/opt 



h2r 2 



-2r x 



N-2 



-2n 



(37) 



In the limit of infinite squeezing, ri,r 2 — > oo, the above 
correlations correspond to a simultaneous eigenstate of 
the relative positions and the total momentum such as 
the cv GHZ states in Eq. @. 

Let us now examine how to experimentally verify the 
genuine multipartite entanglement of the cv GHZ-type 
states (in any case, it may be verified in an operational 
way by doing quantum teleportation between every pair 
of parties with the help of the remaining party Due 
to experimental imperfections, we may assume that the 
entanglement of slightly degraded approximate versions 
of the states generated according to a scheme as in Figs. |l| 
and U is to be verified. We start again with only three 
parties and modes. For a simple check, look at the fol- 




FIG. 2: Verification of genuine tripartite cv entanglement. 
p measurements: directly detecting the p quadratures of all 
three modes and electronically combining them in an ap- 
propriate way. The three-mode tripartite entangled state of 
modes 1,2, and 3 in this figure is produced with three squeez- 
ers and two beam splitters. The parameters gi are the "gains" 
from the conditions in Eq. ( |3^ ) which can be chosen optimally 
(see the text later) . 



lowing set of inequalities, 

I. ([A(ii - x 2 )] 2 ) 

II. ([A(x 2 - x 3 )} 2 ) 



([A(pi +P2 + 93P3)} 2 ) 

([A( 9l pi + p 2 + p 3 )) 2 ) 



> 1 , 

> 1 , 



III. ([A( Xl - x 3 )} 2 ) + ([A(p x + g 2 p 2 + p 3 )} 2 ) > 1 . 

(38) 

On the l.h.s. of condition I., we have hi = —h 2 = gi = 
g 2 = 1 and h 3 — 0, and hence the boundary for the 
total variance in Eq. ([l3|) becomes one with Eq. (JlJ) and 
Eq. (|l6|), but zero with Eq. (14). Similarly, using the 
l.h.s. of condition II., where h 2 = —h 3 = g 2 = <?3 = 1 and 



hi = 0, the boundary is one for Eq. Q and Eq. (|l|), but 
zero for Eq. Jl6|). Finally, the l.h.s of condition III. with 
hi = —h-j = gi = 93 = 1 and h 2 — corresponds to a 
boundary of one in Eq. ( |l4| ) and Eq. ( |l6| ) , and a boundary 
of zero in Eq. (|l5|). Thus, the following statements for 
(at least partially) separable states hold, 



i 

P = Pi,13 ® Pi.2 

i 

P = Z-jVi Pi.2-3 ® Pi.l 



II. and III. , 
I. and II. , 

I. and III. (39) 



The conditions in Eq. (|38|) are necessary for different 
kinds of partial separability. As a result, the violation 
of any pair of inequalities in Eq. (|3^) is sufficient for 
genuine three-party three-mode entanglement. Violat- 
ing only one condition in Eq. ([38|) (for example, con- 
dition I.) means that the total density operator cannot 
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be written in two of the three forms in Eq. (39) (for 
example, neither in the form p — Vi Pi,is ® Pi,2 nor 
in the form p — J^i Vi Pi, 23 ® Pi,i). Using the classifi- 
cation of Ref. |30| , the classes 3 [two- mode biseparable 
states expressible in two of the three forms in Eq. 
4 [three-mode biseparable states expressible in all of the 
three forms in Eq. (|39[) 1, and 5 [fully separable states de- 
scribable by Eq. (\l7\)] are then ruled out. The forms 
of the classes 1 (fully inseparable states) and 2 [one- 
mode biseparable states expressible in one of the three 
forms in Eq. ([59])] remain both possible. In our exam- 
ple with the violation of I., the state might be genuinely 
tripartite entangled or of the partially separable form 
p = J2 i rji pi.yi ® Pi, 3. Eventually, the violation of a sec- 
ond inequality in Eq. (|38l ) (for instance, condition II.) 
negates also the only remaining partially separable form 

(e.g., p = J2iViPi,i2 ®£i 3 )' tnus P rovm g tne Ml insep- 
arability of the state [E0[. Note that even though pure 
and totally symmetric multi-party entangled states are 
always genuinely multipartite entangled |t8] |l9), asym- 
metric pure or mixed entangled three-mode states (e.g., 
from class 2 in Ref. j3(j], the product state of a bipar- 
tite entangled two-mode squeezed state and a vacuum 
state) and symmetric mixed entangled three-mode states 
(like the example for the three-mode biseparable class, 
class 4, given in Ref. [p0| ) do not automatically exhibit 
genuine tripartite entanglement. Due to the violation of 
two conditions in Eq. (|3q), the two loopholes of partial 
separability, mixedness and/or asymmetry, are ruled out. 

The criteria here are only sufficient for full inseparabil- 
ity and hence genuinely tripartite entangled states may 
also satisfy all the conditions in Eq. ( |38| ) (an example will 
be mentioned later). On the other hand, note that we did 
not use the assumption of Gaussian states. The deriva- 
tion of the conditions relies only on the Cauchy-Schwarz 
inequality and Heisenberg's (sum) uncertainty relation. 

Alternatively, one could simply check the known bi- 
partite separability conditions [p5| for pairs of modes, 
i.e., gi = g 2 = 53 = in Eq. (|38D (or using products of 
variances j$8| instead of sums). Again, the statements 
in Eq. ( p9| ) hold. Hence two violations again verify gen- 
uine tripartite entanglement. However, the significance 
of the more general conditions in Eq. (^8|) compared to 
those with g\ — g 2 = g 3 = is that for the cv GHZ-type 
states, as discussed later, the former can be always vio- 
lated for any degree of multi-party entanglement and the 
violations can steadily grow from small towards "perfect" 
(that is all variances of the combinations zero) as the 
three-mode entanglement increases. In contrast, the bi- 
partite conditions with g\ — g 2 — g 3 = may be violated 
for bad three-mode entanglement (small squeezing) and 
satisfied for larger squeezing, thus not always verifying 
genuine tripartite entanglement, and in particular never 
verifying good genuine tripartite entanglement. More- 
over, they might be always violated, but the violations do 
not attain a significant amount (e.g-, three-mode states 
made from one squeezed state |l^, ^9|). Similarly, us- 
ing products of variances pj| instead of sums in Eq. fl3S|) 



with gi = g 2 = .93 = 0, violations may always occur, 
but also only to a certain extent |L8| |l9). In Figs. [l] and 

0, it is shown how to apply the tripartite entanglement 
criteria experimentally using homodyne detectors. 

Let us also discuss the conditions for the N — 4 case 
in more detail. We consider a set of six inequalities, 

1. <[A(i 1 -i 2 )] 2 ) + ([A(p 1 +p 2 + 53 P3+.g4P4)] 2 ) > 1, 

II. ([A(x 2 - x 3 )} 2 ) + ([A( fl ipi + P2 + h + giPi)?) > 1, 

III. ([A{ Xl - x 3 )} 2 ) + ([A{p 1 +g 2 p 2 + p 3 + g4P4)} 2 ) > 1, 
W.([A(x 3 ~x 4 )} 2 ) + ([A(g 1 p 1 +g 2 p 2 +p 3 +p i )} 2 ) > 1, 

V. ([A(x 2 - x 4 )} 2 ) + {[A(gm + p 2 + g 3P3 + p 4 )} 2 ) > 1, 

VI. ([A(ii - i 4 )] 2 > + ([A(pi + g 2 p 2 + g 3 p 3 +Pi)] 2 ) > 1. 

(40) 

The position and momentum variables £1 and pi are 
the quadratures of four electromagnetic modes this time. 
The gi are again arbitrary real parameters. Now the fol- 
lowing statements for (at least partially) separable states 
hold, 



P 



y^ViPi,i23 ® pi,4 -> IV., V., and VI. 



, Vi Pi, 124 <X> Pi, 3 



m Pi,134 <£> Pi,2 



II., III., and IV. 



I.,IL, and V. 



and, 



p = },r)iPi,234 ® Pi,i -> I. ,111., and VI., 

i 

p = ^77iPi,i2 <B> Pi, 34 -> II-, HI-, V., and VI., 

i 

p = y^^r)iPi <13 <8>p»,24 -> I-, II-, IV., and VI., 

i 

p = } j rnpi,i4®Pi,23 -> I., III., IV., and V., 



(41) 



(42) 



and finally, 
P 

P 

P 

P 

P 

P 



m pi.\ 2 d 

i 


§ Pi, 3 $ 


5p i>4 - 


-> all except I., 


^2 Pi A3 e 

i 


5 Pi, 2 i 


5p,,4 - 


-» all except III., 


^ Vi Pi, 14 € 

i 


3 Pi, 2 <i 


5Pi,3 - 


-* all except VI., 


i 


5 Pi,l <i 


§ Pi,4 - 


-> all except II., 


m Pi ,24 <i 

i 


5 pi,i <i 


5^,3 - 


-> all except V., 


X! Vi Pi,34 <5 


5 Pi,i € 


5 Pi,2 ~ 


-* all except IV. 



(43) 
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Note that again the fully separable state, p = JV Vi Pi,i ® 
Pi,2 <S> Pi,3 <8> Pi,4, is included. The above statements can 
be easily confirmed using Eq. (^) for states of the gen- 
eral form Eq. (Pq). The different forms here are p = 
Si % /5i,fcim <S> including p = £\ % ® p i)m (g) 
and p = X)i % Pi.fem ® ^t.Jfii including p = J2iVi Pi,km ® 
Pi,; ® Pi,nj with the two modes m and n always be- 
ing separable. For the combinations u = i m — x n and 
^ = SfcPft +gipi +Pm +Pn, the boundary of the total vari- 
ance is always one. The statements Eq. (flT|), Eq. (42), 
and Eq. (^) become obvious then by considering all pos- 
sible pairs of modes (to, n) of the four modes (fc, I, to, n). 
Note that always when the two modes (to, n) are po- 
tentially entangled, the boundary for the total variance 
drops to zero. 

What kind of violations of the six inequalities in 
Eq. (^) are now sufficient to verify the full inseparabil- 
ity of a four-mode four-party state? The violations must 
rule out any of the partially separable forms in Eq. (41), 
Eq. p^), and Eq. (fl3|). Let us, for example, consider 
violations of the inequalities IV. and V. These viola- 



tions mean that all partially separable forms in Eq. (41), 
Eq. and Eq. (43) are excluded except for the form 

p = Y,iVi Pi,234®Pi,i in Eq. fl4l|). In order to negate this 
form as well a further violation is needed. According to 
Eq. (fill), one of the inequalities I., III., or VI. should be 
violated in addition. Here it is important to realize that 
the conditions IV. and V. do not involve the x quadrature 
of mode 1, but that of all the other modes. The addi- 
tional test via any one of the conditions I., III., or VI., of 
which all contain both quadratures of mode 1, eventually 
provides the missing information about mode 1. Hence 
we learn that three conditions are sufficient here to ver- 
ify the full inseparability of a four-mode four-party state. 
We may choose 

([A(ai - x 2 )} 2 ) + ([A(pi +p2 + g 3 P3 + 9iPi)] 2 ) 
([A(x 2 - x 3 )} 2 ) + {[A(g 1 p 1 +P2+P3+ 54P4)] 2 ) 
{[A(x 3 - x 4 )} 2 ) + ([A(g 1 p 1 + g 2 p 2 + p 3 + p 4 )} 2 ) 



< 1, 

< 1, 

< 1, 



(44) 

which involve both quadratures x and p of all four modes. 
Note that apart from the coefficients gi, these four com- 
binations correspond to those observables measured in a 
four-party cv GHZ state analyzer. Correspondingly, for 
N parties and modes, we may choose the following N — 1 
conditions in terms of effectively N combinations (those 
of an TV-party iV-mode cv GHZ state analyzer), 

([A(xi - x 2 )] 2 ) 

+([A(pi +p 2 + g 3 f>3 H h gNPN)} 2 ) < l, 

([A(x 2 -x 3 )] 2 ) 

+([A(g 1 p 1 +p 2 +P3+ giPi H \- gNPN)} 2 ) < l, 



([A(ijv-i - £n)} 2 ) 
+([A(g 1 p 1 + g 2 p 2 - 




0.2 0.40.60.8 1 

squeezing r 

FIG. 3: Plot of the left-hand-side (total variance) of the con- 
ditions in Eq. (^) for different iV-mode states with quadra- 
ture correlations given by Eq. ([5(i]) and different numbers of 
parties N = 3 and N = 30. The states are those produced 
with one squeezed state (dotted lines with V2 — and r = n) , 
those made from N equally squeezed states (dashed lines with 
r = n = V2), and the unbiased minimum-energy states with 
squeezing ri and r = T2 related as in Eq. (fj4|). 

These conditions are sufficient to verify the full insepa- 
rability (genuine iV-party entanglement) of an iV-party 
iV-mode state. For arbitrary N, the proof relies on the 
fact that in any partially separable form we may always 
select a distinct pair of modes (to, n) which are separable 
in the states i of the convex sum of the density operator. 
Only exploiting that modes to and n are separable, the 
combinations 



N-2 

E 

3=1 



gkjPkj +Pm+Pn , 



(46) 



+ gN-2PN-2 + PN-1 + Pn)} 2 ) 

< 1. (45) 



always yield a boundary of one for the total variance 
using Eq. ( p9| ) for states of the general form Eq. (p6|). By 
taking the pairs of modes (1, 2), (2, 3), (N — 1, N) for 
(m,n), all partially separable forms of the total density 
operator are covered (as demonstrated explicitly for N — 
4) and hence the N— 1 conditions in Eq. ( [45| ) are sufficient 
for genuine TV-party iV-mode inseparability. 

The left-hand-side of the inequalities in Eq. (|4^) is 
shown in Fig.|| for various cv GHZ-type TV-mode states 
differing in the relation between the squeezing rq and 
r 2 [Eq. ©, Eq. (§|), and Eq. ©]. Due to the total 
symmetry of all these states, the left-hand-side of the 
conditions in Eq. (ptq) becomes equal for all conditions 
(assuming gj = g^^). Hence values below the boundary 
1 here mean all inequalities in Eq. ( fl5| ) are satisfied, thus 
indicating genuine iV-party entanglement. In all cases in 
Fig. ||, the optimal coefficients gj = g^ from Eq . (|37|) 
are used to minimize the total variances of Eq. (p6|). If 
N = 30, only for the unbiased states, the conditions are 
always met (for any nonzero squeezing r > 0) and the 
total variances tend to zero for large squeezing. More- 
over, for the same squeezing r, the unbiased states with 
N — 30 drop below the boundary 1 to a greater extent 
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than their unbiased tripartite counterparts. In contrast, 
for the biased states (those with only one squeezer, r 2 = 
and r = r±, and those with N equally squeezed states, 
r = r% = T2), the total variances approach or even ex- 
ceed the boundary 1 as the number of parties grows. The 
example of the states with N equal squeezers also demon- 
strates that there are Gaussian states which are indeed 
genuinely iV-party entangled, but do not satisfy any of 
the conditions in Eq. (|45|). It can be shown, however, tak- 
ing into account the symmetry and purity of the whole 
family of iV-mode states (including those with N equal 
squeezers) that all these states are genuinely multi-party 
entangled for any nonzero squeezing [Jl8|, |l9| . 

Finally, we emphasize that one may use other condi- 
tions too for verifying the genuine multipartite entangle- 
ment of the cv GHZ-type states. Even a single condition 
might be again sufficient. For example, consider the com- 
binations u — 2x\ — {x2 + £3) and v = pi + P2 + P3 for 
three modes. We have [u, v] = 0, and indeed the GHZ- 
type three-mode state becomes a simultaneous eigenstate 
of u and v in the limit of infinite squeezing, r\,r2 —* 00. 
The boundaries of the total variance for these combina- 
tions take on the value one when p — JX Vi Pi-Xi ® Pi, 3 or 
p = J2i Vi Pi. 13 ® Pi. 2 j an d the value two (corresponding 
to the fully separable state) when p = JX rji pi^z ® Pi,i- 
Hence ((Aw) 2 ) p + ((Av) 2 ) p < 1 is sufficient for genuine 
tripartite entanglement. The number of measurements 
required, however, remain the same as for the criteria 
above expressed by N — 1 conditions. In any case, both 
quadratures of all modes must be detected and combined 
in an appropriate way. 

IX. CONCLUSIONS 

In summary, we proposed experimental criteria to de- 
tect genuine multipartite continuous-variable entangle- 
ment. These are expressed in terms of the variances of 
particular combinations of all the quadratures involved. 
The combinations are measurable with only a few simple 
homodyne detections. For Gaussian states, it is then not 
necessary to determine the entire correlation matrix in 



order to confirm the genuine multipartite entanglement. 
Furthermore, the conditions here do not rely on the as- 
sumption of Gaussian states. An experimental confirma- 
tion of the Gaussian character of the state in question 
is therefore not needed cither. Finally, we examined the 
applicability of the conditions to a particular GHZ-type 
class of genuinely multi-party entangled states. These 
states are of Gaussian form, they are totally symmetric 
under exchange of modes, and they have zero cross corre- 
lations between the x and the p quadratures. If they are 
in addition unbiased between the x and the p quadra- 
tures, they always (for any nonzero entanglement) sat- 
isfy the conditions in terms of appropriately chosen lin- 
ear combinations. In the limit of perfect entanglement, 
the variances of the combinations tend to zero for the 
unbiased states and the conditions are perfectly met. 

In an experiment, one normally has approximate a pri- 
ori knowledge about the state to be analyzed. According 
to this a priori knowledge, one can then choose appro- 
priate linear combinations to be measured. It would be 
desirable to know whether there is always, for any given 
multi-party multi-mode state, a single optimal condition 
to verify its genuine multipartite entanglement and how 
to constructively derive this condition. Inferring from 
the results here, such a condition may always exist and 
the corresponding linear combinations must contain both 
quadratures of all modes with optimized coefficients hi 
and gi. A possible approach to this question is in terms 
of so-called entanglement witnesses [^6[ One may 

then interpret the inequalities for the total variances 
as quantum expectation values of Hermitian operators 
which take on negative values when they witness some 
kind of partial inseparability. 
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